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Abstract 

The stationary phase method is often employed for computing tunnehng phase times of 
analytically-continuous gaussian or infinite-bandwidth step pulses which collide with a potential 
barrier. The indiscriminate utilization of this method without considering the barrier boundary 
effects leads to some misconceptions in the interpretation of the phase times. After reexamining 
the above barrier diffusion problem where we notice the wave packet collision necessarily leads to 
the possibility of multiple reflected and transmitted wave packets, we study the phase times for 
tunneling/reflecting particles in a framework where an idea of multiple wave packet decomposition 
is recovered. To partially overcome the analytical incongruities which rise up when tunneling phase 
time expressions are obtained, we present a theoretical exercise involving a symmetrical collision 
between two identical wave packets and a one dimensional squared potential barrier where the 
scattered wave packets can be recomposed by summing the amplitudes of simultaneously reflected 
and transmitted waves. 
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I. INTRODUCTION 



The analytical methods utilized for reproducing the wave packet collision with a po- 



tential barrier have been wide 
tunneling phenomena 



y discussed in the context of (one dimensional) scattering and 
3|]. In this extensively explored scenario, the stationary phase 
method (SPM) is the simplest and the most common approximation employed for obtaining 
the group velocity of a wave packet in a quantum scattering process which represents the 
collision of a particle with a square (rectangular) barrier of potential. First introduced to 
physics by Stokes and Kelvin jj], the SPM essentially provides an approximate value for 
the maximum of an integral in an analytical procedure which enables us to understand sev- 
eral subtleties of interesting quantum phenomena, such as tunneling [l, [s], Q], resonances 



iartman 



ll|. It is 



[3], incidence- reflection and incidence-transmission interferences 
effect [8| and the correspondent superluminal transition time interpretation 

an 

also often quoted in the context of testing different theories for temporal quantities such as 
arrival, dwell and delay times y, [l2] and the asymptotic behavior at long times [ill, isj ]. 

The study presented in this manuscript concerns with some limitations on the use of the 
SPM for obtaining the right position of a propagating wave packet subject to a very simple 
and common potential configuration. The limitations are not concerned with the domain 
of the method, i. e. with the mathematical basis and the analytical applicability of the 
method, they focus on the indiscriminate application of the method to time-dependent scat- 
tering problems and, in particular, for deriving phase times. The time- dependent scattering 
of a wave-packet by a potential barrier in the linear Schroedinger equation is interesting 
because it reconciles two types of wave theories. Those interested in linear theory often look 
at stationary (frequency-domain) solutions of a single frequency, whereas those interested in 
nonlinear theory usually look at time-dependent pulses. To implement the accurate analysis 
with the SPM, we work with the hypothesis of adding up the Fourier modes from the station- 
ary solution to observe the evolution of a linear wave packet on collision with an obstacle. In 
the first part of the manuscript, the analysis of the applicability of the SPM is concentrated 
on the investigation of the non-relativistic diffusion of an incoming single wave packet with 
energy spectrum totally localized above the potential barrier. We are particularly interested 
in obtaining an analytical description of the above barrier diffusion problem by pointing 
out some interpretive problems which concern with non-conservation of probability during 
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the collision process. In the second part of our analysis, the correct quantification of the 
analytical incongruities which restrict the applicability of the principle of stationary phase 
is extended to the study a tunneling particle where, eventually, the superluminal interpreta- 
tion of transition times can be ruined. In order to organize the ideas which are discussed in 
the manuscript, we have stated the following division. In section II we discuss the general 
applicability of the SPM in a wave packet propagation problem by introducing an analyt- 
ical description of the above barrier diffusion problem where we reobtain the solution for 
the apparent physical inconsistency which concerns with the non-conservation of probab ility 
during the collision process. By recurring to a wave packet multiple peak decomposition l4] , 
we can decompose the original colliding wave packet into multiple reflected and transmitted 
components by means of which the conservation of probabilities is recovered during the scat- 
tering process. The analytic expressions for all the reflected, transmitted and intermediary 
components are obtained and the validity of them is discussed in terms of an illustrative 
example where the analytic approximation and the reflects exactly the numerical results. 
The physical framework where the study of a tunneling process is embedded is permeated by 
theoretical constructions presented in section III. It involves analytically-continuous gaus- 
sian pulses which due to its not well-deflned front lead to ambiguities in the interpretation 
of the wave packet speed of propagation. Beside it, inflnite bandwidth signals cannot prop- 
agate through any real physical medium (whose transfer function is therefore finite) without 
pulse distortion, which also leads to ambiguities in determining propagation velocity during 
the tunneling process. Moreover, some of the barrier transposing time definitions lead, in 
tunneling time conditions, to very short times which can even become negative where they 
may seem to contradict simple concepts of causality. To partially overcome these incom- 
patibilities, in section IV, we discuss a theoretical exercise involving a symmetrical collision 
with a one dimensional square potential where we recompose the scattered wave packet by 
summing the amplitudes of the reflected and transmitted waves in the scope of what we 
classify as a multiple peak decomposition analysis [l^ where the conditions for the SPM 
applicability are totally recovered. Finally, we draw our conclusions in section V. 
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II. THE ABOVE BARRIER DIFFUSION PROBLEM 



The principle of stationary phase to describe the group velocity is well-established, and 
was employed by Sommerfeld and Brillouin in their early studies of wave propagation and 
group velocit;^ 15|]. The method has in time become a standard tool for several theoretical 
applications [16] adopted not only by physicists but also by biologists, economists etc. It 
can also be used in a statistical sense, such that the most likely events tend to be associated 
with slowly- varying phase variation in the frequency domain, and unlikely events tend to be 
associated with rapidly varying phase with frequency. 

To illustrate the applicability of the SPM into a physical problem it should be instructive 
to assume an equivalence between the complex variable z and a real one k so that an integral 
which represents a free quantum particle propagation can be identified by the wave packet 
solution of the unidimensional Schroedinger equation. 



.. 9 , , , 

with V{x) = 0, described in terms of the integral 

dk 



+ V{x) 



G{k, ko) exp [—1 {E{k) t + kx — k Xq)], 



(1) 



(2) 



where we have set ^ = 1 and the dispersion relation as E{k) = A;^/(2m). The integral ([2]) 
can therefore be estimated by finding the place where the phase has a vanishing derivative, 
evaluating (approximately) the integral in the neighborhood of this point. The movement of 
the peak coordinate of the wave packet ip{x^t) can be obtained by imposing the stationary 
phase condition 



^[Et- k{x-Xo)] 



k=ko 



Xq -\ t, 

m 



(3) 



which means that the maximum of the wave packet propagates with a velocity v = In 
fact, we ratify this result when we explicitly calculate the integral of Eq. ([2]) by introducing 
a gaussian momentum distribution given by 

a^{k — ko) 



G{k,ko) = g{k - ko) = exp 
In this case, the result of the integration ([2]) give us 



(4) 



ll}{x,t) = Lp[x - Xo,t], 



(5) 
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where 



na 
2 



1 + 



if 



exp 



fco i 
m 



(1 + ^) 



arctan {j^) + i(A;o2; - Eot) 



(6) 



which evidently ratifies the result of Eq. ([3]). 

Meanwhile, the method leads to some new interpretive discussion when the momen- 
tum distribution becomes a complex function, i. e. when G{k, kg) has to be written as 
\G{k, ko)\ exp [i X[k)]. In this case, the stationary phase condition for a free wave packet 
represented by Eq. ([3]) will be modified by the presence of an additional phase \{k). By 
assuming that X{k) can be expanded in the neighborhood of kg, i. e. 

dX{k) 



X{k) « A(fco) + {k- fco) 

the new stationary phase condition should give 

dX{k) 



dk 



k=ko 



dk 



k=ko 



+ —t. 

m 



(7) 



(8) 



As it is commonly presented in several textbooks of quantum mechanics 17|, one can be 
persuaded to take the stationary phase condition as a necessary and sufficient statement 
for applying the SPM. Nevertheless, although the principle of stationary phase can be used 
to predict the most likely outcome of an experiment, it does not exclude the possibility of 
alternative outcomes. In fact, the results of the SPM depend critically upon the manipulation 
of the generic amplitude G{k, kg) prior to the application of the method. To clear up this 
point, let us firstly analyze the diffusion of a wave packet by a squared potential K when 
the incoming energies are higher than Vq. The stationary wave solution of the Schroedinger 
equation ([T]) computed by considering the potential barrier described by 

if a: < region I, 

Vo if < X < L region II, (9) 
if X > L region III, 

can be decomposed into different wave functions for each interval of x, i. e. 

^{k,x) = (t)i{k,x) + (t)ii{k,x) + (t)iji{k,x), 

(pinc{k,x) + (t)R{k,x) (j)^{k,x) + (t>p{k,x) (t>Tik,x) 



(10) 
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where 



4)inc{k,x) = exp[ikx], 

(puikjx) = R{k)eyip[—ikx], 

4)a{k,x) = a(A;) exp [i g x], 

(j)i3{k,x) = I3{k)exp[—iqx], 

(prikjx) = T{k) eyip[i k x], (11) 

with q = (k^ — w'^)^ and w = a/2 m V^. By solving the constraint equations at x = and 
X = L we obtain 

'k{k + qy 



a{k) 

m 

R{k) 
T{k) 



Hk) 
'k{k-q) 



Hk) 



exp [i Q{k) — iqL], 
exp [i Q(k) + iq L], 
sin [q L] exp [iQ{k)], 



2kq 
Hk) 



exp [iQ{k) — ik L], 



where 



^(/fc) = {Ae cos" [q L] + {k^ + qysin^[q L\Y 



and 



Q{k) 



arctan 



k"^ + q 
2kq 



■ tan [q L] 



(12) 



(13) 



(14) 



The exphcit expression for the correspondent propagating wave packets can be obtained by 
solving the integrals like 

dk 



tpf{x,t) 



g{k — ko) 4>f{k, x) exp [—lEt], 



(15) 



with / = g, /3, i?, T. As first approximation which is commonly used in quantum mechanics 



textbooks 



17l |. we obtain the analytical formulas to these integrals by assuming the momen- 



tum distribution g{k — k^) is sufficiently sharply peaked around a maximum point > w. 
In this case, the integration can be extended from [w, oo] to [— oo, oo] without modifying 
the final result. However, in the sense we are investigating, as we shall demonstrate in the 
following, it is erroneously assumed that the /c-dependent phase terms Q{k) and q can be 
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approximately described by a series expansion around k = np to the first order term, i. 
e. 



e{k)^e{ko) + {k-k,)Q'{ko) and q ^ + {k - ko) 



dk 



(16) 



k=ko 



where 



{k^ + q^)k^qL- {k^ - q^y sin [q L] cos [q L] 



dq 

Tk 



ko 



k=ko 



Ak^q^ + (A;2 - q^fsin^ [qL] 
At the same time, when k is approximated by ko, it is assumed that the modulating amplitude 
\(pf{k,x) \ can be put out of the integral giving the following results. 



1plncix,t) ip[x - Xo,t], 



dk 



g{k — ko) exp [—iEt — i k{x + Xq) + +i{k — /c„)0'(A;o)] 



i'^{x,t) PS a{ko) 



R{ko) exp [-1 ko Q'iko)] (f [-x - Xq + @'{ko),t], 
dk 



g{k - ko) exp \-i {Et + kxo-qox)+i{k-ko)(^^{x-L)+ &{ko] 



Q(A;o)exp iqoX - i-^ {x - L) - ikoQ'{ko) if -^{x - L) - Xo + &' {ko),t 



ipp{x,t) K P{ko) 



dk 



g{k-ko)exp -i{Et + kxo + qox)-i{k-ko)(^ix-L)-Q'{ko] 



l3{ko) exp -iqoX + i-^{x- L)-i h0{ko] 



^{x-L)-xo + e'{ko),t 



<?0 



'ipTix,t) T{ko) 



dk 



g{k — ko) exp [—iEt + i k{x — Xo) + i{k — ko) {Q'{ko) — L)] 



T{ko) exp [-1 koQ'{ko)] (p[x - Xo - L + Q'{ko),t]. 



(17) 



By observing the gaussian shape of given by Eq. ([6]), one can easily (but wrongly!) 

identify the position of the peak of the wave packets. The times corresponding to the position 
X of the incident and reflected wave packet peaks would be respectively given by 

x + Xo- @'{k)' 



tlncyX) 



where Vk — '^^ — ^ 



Vk 



and tfl:(x) 



k=ko 



Vk 



(18) 



k=ko 



dk ~ m- Only X < is physical in this result since these waves, by definition, lie 
in region I. Since the phase of the incoming wave contains only the plane wave factors, i.e. it 
is devoid of Q{k), the incoming peak reaches the barrier at x = at time t/„c(0) = — (a^o/^fco) 
(neglecting interference effects). The presence of the phase term Q'{ko) = {dQ{k) / dk)\k=k^ 
for tfl,(x) would imply a time delay, or a time advance, which depends on the signal of G'(/co), 
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for the reflected wave with respect to the incident one, i.e. 

e'(A;) 



tB.{0) = tinciO) 



Vk 



(19) 



k=ko 



t is analogous to what happens for the step potential tunneling penetration when E < Vq 
~r\ . When < a; < L in region II, the times for the "intermediary" a and (3 wave packet 



peaks would be given by 

'{x-L) {x,-Q'{k)) 



Vk 



and ti3{x) 



k=kn 



(x-L) 



jxo - e'jk)) 

Vk 



(20) 



And finally, when x > L in region III, the peak of the transmitted wave packet would be 
written as 



trix) 



X- Xo- L + e'(fc) 



Vk 



(21) 



k=ko 



The above results fll8H21l) could be directly obtained by applying the SPM to the wave 
functions ipf{x,t) expressed in Eq. ( fTSl) . where the stationary phases for each decomposed 
wave component / would be given by 

'&lncix,t,k) = -Et + k{x - Xo), 

'&R{x,t,k) = -Et - k{x + Xo) + Q{k), 

'd^{x,t,k) = -Et - kxo + q{x - L) + Q{k), 

'df3{x,t,k) = -Et - kxo- q{x - L) + Q{k), 

^T{x,t,k) = -Et + k{x-Xo-L) + @{k). (22) 

Meanwhile, these time-dependencies must be carefully interpreted. In the Fig. ([2]), we il- 
lustrate the wave packet diffusion described analytically by the results of Eq. (fT7|) . The 
"pictures" display the wave function in the proximity of the barrier for suitably chosen 
times. By taking separately the right (a) and left {(3) moving components in region II, 
independent of the value of Q{k), we can observe that 

L 



At, = t^{L) - t^{0) = At I, = t^(0) - t^iL) 



(23) 



k=ko 



which should corresponds to transit time values "classically" expected. However, by con- 
sidering the information carried by the wave packet peaks, we would have a complete time 
discontinuity at x = represented by Eq. f|T9l) and by the fact that 

L Q'{k)' 



Vk Vq 



Vk 



+ tluM, 



(24) 



k=ko 



i. e. the a wave could appear in region II at a time tc«(0) before (or after) the incident wave 
having arrived at x = 0. Both results expressed by Eqs. f|T9|) and flMl) which are illustrated 
in Fig.[2]should imply the non-conservation of probabilities. It is, by principle, unacceptable. 

To be more accurate with this analysis and clear up some dubious points, it should be 
more convenient to verify two simple particular cases. In order to simplify the calculations, 
we set Xq = and we choose go L = ravr in = 1, 2, 3, ...) so that we can write. 



.F(fco) = 2fcogo, e(A;o) = 0, ii(fco) = 0, T{k,) = \ (25) 

and 

9>'ih\ = - 

so that 



e'(fc) = ^^1±A > ^ (26) 



tM = ^ " ' > 0. (27) 

Otherwise, if we choose q^L = {n + |)7r, we shall obtain, 

Hko) = K + <ll m = I, R{ko) = T(fco) = 1^ exp \i (I - L)h8) 

and 



Qo K + ql qo 

which gives 



e'{k) = H^p^ < UlL (29) 



\2 



, , mL (kn — Qo) 

= ^ ° J'^ < 0. 30 

qo kl + ql 

The latter negative value corroborates with the absurd possibility of the peak associated 
with the a wave appears before the peak of the incident wave packet arrival at a; = 0. In 
fact, we must observe that if we proceed with a more general analysis, the phase derivative 
&{k) does not present a regular behavior. It oscillates very rapidly, as we can observe in 
Fig.([T]) so that the SPM cannot he applied! 

Other inc ong ruities in the naive application of the SPM presented above have been 
pointed out 14|. Differently from the analytical analysis presented above, the numerical 
simulations of a wave packet diffusion above the potential barrier shows the appearance 
of multiple peaks due to the two reflection points at x = and x = L (see the figures). 
In fact, numerical calculations automatically conserve probabilities, at least to within the 
numerical errors. This observation suggest a new analysis and subsequent interpretation of 
the ambiguities presented in the previous section. In order to correctly apply the SPM and 



accurately obtain the position of the peak of the propagating wave packet with an accurate 
time dependence, we are constrict to solve the continuity constraints of the Schroedinger 
equation at each potential discontinuity point x = and x = L hj considering multiple 
successive reflections and transmissions. The phenomenon was already described in and 
consists in an incoming wave of unitary amplitude with momentum distribution centered 
at fco which reaches the interface a; = where, at an instant t = t^, is decomposed into a 
reflected wave component of amplitude (see the diagram) and a transmitted wave com- 
ponent of amplitude Oi. The transmitted wave continues propagating until it reaches x = L 
where, at t = ti, it is now decomposed into a reflected wave component of amplitude /3i 
and a transmitted wave component of amplitude T^. That new reflected wave component 
comes back to a; = where, at t = ta, it is again decomposed into another reflected wave 
component of amplitude and transmitted wave component of amplitude R2 which will be 
added to R^. This iterative process continues for infinity times as we can see in the following 
diagram. 



exp[i kx] + R-i exp[— i k x] 
R2 exp[— i k x] 



«! exp[z q x\ 

mexp[igx] + /3iexp[— zgx] 
: + : 



aiex.p[iqx] + /3iexp[— zgx] 
a2exp[iqx] + /32exp[— igx] 



+ 



R„exp[—ikx] a„ exp[i gx] + /3„_i exp[— i q x] a„exp[iqx] + P„exp[—iqx 
X = x = L 

The continuity constraints over ipi^, t) for each potential step at a; = and x = L determine 
the coefficients 

k-q 2 k „ 2k{q-k) 



Ti exp[i k x] 
T2 exp[i k x] 

T„ exp[i k x] 



k + q' 
Akq 



k + q' 
exp[i {q — k) L] 



/3i 



R, 



(k + q) 

and establish the recurrence relations 



(k + qy 

— a\ pi, 



exp [2 i g L ] , 



R„ 



R 



n + l 



/3n 



T 



k 



k + q 



exp[2 iqL] 



n 



1,2, 



(31) 



(32) 
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which allow us to write the sum of the coefficients a„, and T„ as 




71=1 '- ^ ' 

OO r / ) \ 2 n — 1 





T = j;r„ =T, l-(^) exp[2zgL] 

rj,= i L \ y / 



(33) 



These sums reproduce exactly the expressions in Eq. ( fT2l) . In this form the interpretation 
is easy. represents the ffist reflected wave (it has no time delay since it is real). R2 
represents the second reflected wave and, as a consequence of the continuity condition at 
X = 0, it is the sum, in region II, of the ffist left-going wave (/3i) and the second right-going 
amplitude (02), i.e. 



This structure is that given by considering two "step functions" back-to-back. Thus at each 
interface the "reflected" and "transmitted" waves are instantaneous i. e. without any delay 
time. 

Now we can calculate the time at which each transmitted or reflected wave appears by 
applying the SPM for each component of the total transmitted T or reflected R coefficients. 
It will give us the recurrence relation 



which coincides with the "classically" predicted value for the velocity of the particle above 
the barrier. Indeed the SPM applied separately to each term in the above series expansion 
for R yields delay multiples of 2 {dq/dE)\^^^^ L = 2(m/go)-^- This agrees perfectly with 
the fact that since the peak momentum in region II is go, the a and (3 waves have group 
velocities of qo/m and hence transit times (one way) of {m/qo) L. The ffist transmitted peak 
appears (according to this version of the SPM) after a time {m/qo) L, in perfect agreement 
with the above interpretation. 

Is this compatible with probability conservation? It is because of the following identity 



R2 = a2 + i3i = - ai f3i . 



00 




(35) 



n=i 
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This result is by no means obvious since it coexists with the well-known result, from the 
plane wave analysis, 

oo oo 

\R\' + \T\' = \J2Rr.\' + \J2T„\' = l. (36) 

n=i n=i 

To conclude, we can state that the conditions for applying the SPM depend upon the correct 
manipulation of the amplitude. A posteriori this seems obvious but the point is that unless 
we know the number of separate peaks the SPM is ambiguous. There is, however, a converse 
to this question. If the modulating function is such that two or more wave packets overlap, 
then we cannot treat them separately without considering all the interference effects. The 
above barrier analysis is simply a particular example of this ambiguity. 

In order to obtain an analytic formula for the wave packet propagation, we consider 
the gaussian momentum distribution given by Eq. (jlj). Using the linear approximation for 
the momentum q given by (fT6|) and considering the expressions in (133|) resulting from the 
multiple peak decomposition in (fT5|) . again, with the same pertinent approximations used 
for obtaining (1171) . we can analytically construct the following simplified expressions 

tpi{x,t) = ^p[x - Xo, t], 



kp-qo 
ko+qo 



^[-x- Xo, t] + 



1pc.{x,t) 

i>f3{x,t) 
i^T{x,t) 



4fcogo(go-fco) 
(fco+go)=* 



qo 



'ip[-x-Xo + 2{n + l)^L, t] 



-i^x 
qo 



E 



kp-qo 
kp+qo 



[- 



exp L 



"1 
qp 



2kp(qp-kp) 



¥■ exp 



(fco+<Jo) 



(kp+qpY ^^P 



kp-qp 



qp 



ip[i2nL + 2L-x)^-Xo, t], 



qo 



^1 E (fcf -p 



qo 



ip[x-Xo-L+i2n + l)^L,t] 



(37) 



which allow us to visualize the multiple peak decomposition illustrated in Fig. [31 By com- 
paring with the exact numerical results, the analytic approximations expressed in flTTI) will 
be valid only under certain restrictive conditions. To explain such a statement, when we 
analytically solve the integrals like (fT5|) . the exponential exp [2iqL] which appears in the 
recurrence relations (l32ll are approximated by 

^0 1 



exp [2i Qo L + 2i {k — ko) 



Qo' 



(38) 
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In this case, the above analytic expressions can be obtained only when the exponential 
function does not oscillate in the interval of relevance of the envelop gaussian function 
g{k — ko), i. e. when AqL < tt, which can be expressed in terms of the wave packet width 

as 



AkL^^-<7r. (39) 

k=ko 9o a 



where we have used A A; ^ 1/a. The above constraint implicitly carries a very peculiar 
character: the multiple peak decomposition can also be evidenced when the wave packet 
width a is larger than the potential barrier width L, i. e. when we assume the relation ko/qo 
is in the interval 1 < fco/o'o < vr we can have propagating wave packets with a > L satisfying 
the requirements for the multiple peak resolution ([3] 



III. LIMITATIONS ON THE TUNNELING TIME ANALYSIS 

On the theoretical front, people have tried to introduce quantities that have the dimension 
of time and can somehow be associated with the passage of the particle through the barrier 



or, strictly speaking, with the definition of the tunneling time. Since a 



efforts have led to the introduction of several time definitions 



24 



0, Q, Q, 



19 



ong time these 



20 



21 



22 



251], some of which are completely unrelated to the others, which can be organized into 

three groups: (1) The first group comprises a time-dependent description in terms of wave 

packets where some features of an incident wave packet and the comparable features of the 

n 

transmitted packet are utilized to describe a delay as tunneling time [1] . (2) In the second 
group the tunneling times are computed with basis on averages over a set of kinematical 
paths, whose distribution is supposed to describe the particle motion inside a barrier, i. 
e. Feynman paths are used like real paths to calculate an average tunneling time with 
the weighting function exp [i S x{t)/h], where S is the action associated with the path x{t) 
(where x(t) represents the Feynman paths initiated from a point on the left of the barrier 



and ending at another 
the Bohm approach 



27 



joint on the right of it [26]). The Wigner distribution paths ^], and 
] are included in this group. (3) In the third group we notice the 
introduction of a new degree of freedom, constituting a physical clock for the measurements 
of tunneling times. Separately, it stands by itself the dwell time approach. The time related 
to the tunneling process is defined by the interval during which the incident flux has to exist 
and act, to provide the expected accumulated particle storage, inside the barrier j^. The 
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methods with a Larmor Cock Q or an oscillating barrier arc comprised by this gxonp. 

There is no general agreement [9|, |25| among the above definitions about the meaning of 
tunnehng times (some of the proposed tunnehng times are actually traversal times, while 
others seem to represent in reality only the spread of their distributions) and about which, if 
any, of them is the proper tunneling time, in particular, due to the following reasons [25] : (a) 
the problem of defining tunneling times is closely connected with the more general definition 
of the quantum-collision duration, and therefore with the fundamental fact that in quantum 
mechanics, time enters as a parameter rather than an observable to which an operator can 
be assigned (b) the motion of particles inside a potential barrier is a quantum phenomenon, 
that till now has been devoid of any direct classical limit; (c) there are essential differences 
among the initial, boundary and external conditions assumed within the various definitions 
proposed in the literature; those differences have not been sufficiently analyzed yet. In 
particular, the study of tunneling mechanisms is embedded by theoretical constructions 
involving analytically-continuous gaussian, or infinite-bandwidth step pulses to examine the 
tunneling process. Nevertheless, such holomorphic functions do not have a well-defined 
front in a manner that the interpretation of the wave packet speed of propagation becomes 
ambiguous. Moreover, infinite bandwidth signals cannot propagate through any real physical 
medium (whose transfer function is therefore finite) without pulse distortion, which also leads 
to ambiguities in determining the propagation velocity during the tunneling process. For 
instance, some of the barrier traversal time definitions lead, in tunneling time conditions, to 
very short times which can even become negative, precipitately inducing to an interpretation 
of violation of simple concepts of causality. Otherwise, negative speeds do not seem to create 
problems with causality, since they were predicted both within special relativity and within 



quantum mechanics 



23j. A possible explanation of the time advancements related to the 



negative speeds can come, in any case, from consideration of the very rapid spreading of the 
initial and transmitted wave packets for large momentum distribution widths. Due to the 
similarities between tunneling (quantum) packets and evanescent (classical) waves, exactly 
the same phenomena are to be expected in the case of classical barriers (we can mention 
the analogy between the stationary Helmholtz equation for an electromagnetic wave packet 
- in a waveguide, for instance - in the presence of a classical barrier and the stationary 
Schroedinger equation, in the presence of a potential barrier . The existence of 

such negative times is predicted by relativity itself, on the basis of the ordinary postulates 
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251], and they appear to have been experimentally detected in many works 3ll. |32|. 



In this extensive 
phase times 



33 



y explored scenario, the first group quoted above contains the so-called 
34 | which are obtained when the stationary phase method (SPM) is em- 
ployed for obtaining the times related to the motion of the wave packet spatial centroid 
which adopts averages over fiuxes pointing in a well-defined direction only, and has recourse 
to a quantum operator for time [25|. Generically speaking, the SPM essentially enables us 
to parameterize some subtleties of several quantum phenomena, such as tunneling |l, ^, Q], 



resonances 



11 



as well as the 
2511 ■ In fact, it 



[31, incidence-refiection and incidence-transmission interferences 
Hartman Effect [8[ and its superluminal traversal time interpretation 
is the simplest and commonest approximation method for describing the group velocity of 
a wave packet in a quantum scattering process represented by a collision of a particle with 



a potential barrier 



14 



19 



25 



In the following, our attention is particularly concentrated on some limitations on the 
use of the SPM for deriving tunneling times for which we furnish an accurate quantification 
of the analytical incongruities which restrict the applicability of this method. We introduce 
a theoretical construction involving a symmetrical collision with a unidimensional square 
potential where the scattered wave packets can be reconstructed by summing the amplitudes 
of the refiected and transmitted waves in the scope of what we denominate a multiple 



peak decomposition analysis 



14| in a manner that the analytical conditions for the SPM 



applicability are totally recovered. 

Generically speaking, the SPM can be successfully utilized for describing the movement 
of the center of a wave packet constructed in terms of a momentum distribution g{k — ko) 
which has a pronounced peak around k^. By assuming the phase which characterizes the 
propagation varies sufficiently smoothly around the maximum of g{k — ko), the stationary 
phase condition enable us to calculate the position of the peak of the wave packet (highest 
probability region to find the propagating particle). With regard to the tunneling effect, the 
method is indiscriminately applied to find the position of a wave packet which traverse a 
potential barrier. 
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For the case we consider the potential barrier 



V{x) 



Vo 



X G [-L/2, L/2] 



(40) 



X ^ [-L/2, L/2] 

it is well known that the transmitted wave packet solution {x > L/2) calculated by means 
of the Schroedinger formalism is given by 17 1 



r dk 

V{x, t) = j i^g{k- ko) \T{k, L)\ exp 



ik(x- L/2) -i t + i Q(k, L) 

2 m 



where, in case of tunneling, the transmitted amplitude is written as 

,,4 



\T{k,L) 



1 + 



w 



4A;2 

and the phase shift is described in terms of 

2k^-w^ 



sinh^ [p{k) L] 



Q{k, L) = arctan 



kp{k) 



tanh [p{k) L] \ , 



(41) 



(42) 



(43) 



for which we have made explicit the dependence on the barrier length L and we have adopted 
p{k) = [w"^ — k"^)'^ with w = (2mV"o)^ and h = 1. Without thinking over an eventual 
distortion that \T{k,L)\ causes to the supposedly symmetric function g{k — ko), when the 
stationary phase condition is applied to the phase of Eq. fH2]) we obtain 



±L{x-L/2)-^t + e{k,L) 



x-L/2 



t + 



d@{k,L) 



0. 



(44) 



m dk 

The above result is frequently adopted for calculating the transit time It of a transmitted 
wave packet when its peak emerges at x = L/2, 



tT 



m 



de{k,a^^^) 



k„ 



dk 



2m L \ sinh (a) cosh (a) — (2 k^ 



w 



2)A;2 a 

J max 



. a 



Akl 
•2 - A; 



^Lx) + sinh (a) 



(45) 



where we have defined the parameter a = {w^ — k^^^J ^ L. The concept of opaque limit is 
introduced when we assume that k^^^ is independent of L and then we make a tends to oo 



In this case, the transit time can be rewritten as 

OL 2m 

= k oik ) ■ 

•^max /-'V max/ 

In the literature, the value of fc^ax is frequently approximated by kg, the maximum of g{k—ko), 
which, in fact, does not depend on L and could lead us to the transmission time superluminal 
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interpretation 



an 



35| of (H6ll . To clear up this point, we notice that when we take the 
opaque hmit with L going to oo and w fixed as well as with w going to oo and L fixed, with 
kg < w in both cases, the expression (H^ suggests times corresponding to a transmission 
process performed with velocities higher than c [11,]. Such a superluminal interpretation 
was extended to the study of quantum tunneling through two successive barriers separated 
by a free region where it was theoretically verified that the total traversal time does not 
depend not only on the barrier widths, but also on the free region extension [oj] . Besides, in 
a subsequent analysis, the same technique was applied to a problem with multiple successive 



barrier where the tunneling process was designated as a highly non-local phenomenon [35 1. 

It would be perfectly acceptable to consider k^^^ = k^ for the application of the stationary 
phase condition when the momentum distribution g{k — ki^) centered at k^ is not modified by 
any bound condition. This is the case of the incident wave packet before colliding with the 
potential barrier. Our criticism concern, however, with the way of obtaining all the above 
results for the transmitted wave packet. It does not have taken into account the bounds and 
enhancements imposed by the analytical form of the transmission coefficient. 

Hartman himself had already noticed that the 

" Discussion of the time dependent behavior of a wave packet is complicated by its dif- 
fuse or spreading nature; however, the position of the peak of a symmetrical packet can be 
described with some precision" 

To perform the correct analysis, we should calculate the right value of fc^ax to be substi- 
tuted in Eq. ( I45l) before taking the opaque limit. We are obliged to consider the relevant 
amplitude for the transmitted wave as the product of a symmetric momentum distribution 
g{k — fco) which describes the incoming wave packet by the modulus of the transmission am- 
plitude T{k, L) which is a crescent function of k. The maximum of this product representing 
the transmission modulating function would be given by the solution of the equation 

'g'ik-k,) , |r(A:,L)|'" 



g{k-ko) \T{k,L)\ 



(47) 



gik-ko) \T{k,L)\ 

Obviously, the peak of the modified momentum distribution is shifted to the right of ko so 
that fc^ax have to be found in the interval J/co;"?^!- Moreover, we confirm that fc^^x presents 
an implicit dependence on L as we can demonstrate by the numerical results presented in 
Table 1 where we have found the maximum of g{k — k^) \T{k, L)\ by assuming g{k — k^) is 
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a gaussian function almost completely comprised in the interval [0,?/;] given by 



g{k - ko) = ( ^ ) exp 



,2 



a" \ " a^{k- kaY 



(48) 



4 

By increasing the value of L with respect to a, the value of k^^^ obtained from the numerical 
calculations to be substituted in Eq. (H5l) also increases until L reaches certain values for 
which the modified momentum distribution becomes unavoidably distorted. In this case, 
the relevant values of k are concentrated around the upper bound value w. We shall show 
in the following that the value of L which starts to distort the momentum distribution can 
be analytically obtained in terms of a. 

Now, if we take the opaque limit of a by fixing L and increasing w, the above results 
immediately ruin the superluminal interpretation upon the result of Eq. (H5l) since tends 
to oo when k^^^ is substituted by w. 

Otherwise, when w is fixed and L tends to oo, the parameter a calculated at k = w 

becomes indeterminate. The transit time tr still tends to oo but now it exhibits a peculiar 

dependence on L which can be easily observed by defining the auxiliary function 

sinh (a) cosh (a) — a 
sinh^ (a) 

which allow us to write 



Gia) = (49) 



t<^='^^G{a). (50) 
w a 

When a S> 1, the transmission time always assume infinite values with an asymptotic 
dependence on {w"^ — k"^ 



.2 U2\-^^ 



^ — ^ OO. (51) 

w {w"^ — k^) 2 

Only when a tends to we have an explicit linear dependence on L given by 

2mL fG{a)\ AmL 



t% = lim <^ -^-^ } = (52) 

w a^o y a ] 

In addition to the above results, the transmitted wave must be carefully studied in terms 
of the rapport between the barrier extension L and the wave packet width a. For very thin 
barriers, i. e. when L is much smaller than a, the modified transmitted wave packet presents 
substantially the same form of the incident one. For thicker barriers, but yet with L < a, 
the peak of the gaussian wave packet modulated by the transmission coefficient is shifted 
to higher energy values, i. e. k^^^ > increases with L. For very thick barriers, i. e. 
when L > a, we are able to observe that the form of the transmitted wave packet is badly 
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distorted with the greatest contribution coming from the Fourier components corresponding 
to the energy w just above the top of the barrier in a kind of filter effect. We observe that 
the quoted distortion starts to appear when the modulated momentum distribution presents 
a local maximal point dX k = w which occurs when 

±[g{k-k,)\T{k,L)\] >0. (53) 

"'^ k=w 

Since the derivative of the gaussian function g{k — k^) is negative aX k = w, the Eq. (153|) 
gives the relation 



g'{w - ko) 
— — < hm 

g[w — ko) k^w 



T'{k,L) 
Tik,L) 



< (54) 

1 + ^' " 



which effectively represents the inequality 



-{w-K)<^ =, L>^-a[l--j. (55) 

Due to the filter effect, the amplitude of the transmitted wave is essentially composed 
by the plane wave components of the front tail of the incoming wave packet which reaches 
the first barrier interface before the peak arrival. Meanwhile, only whether we had cut 
the momentum distribution off at a value of k smaller than w, i. e. k ^ [1 — 6)w, the 
superluminal interpretation of the transition time (H6|) could be recovered. In this case, 
independently of the way as a tends to oo, the value assumed by the transit time would 
be approximated by t^ ^ 2m/w6 which is a finite quantity. Such a finite value would 
confirm the hypothesis of superluminality. However, the cut off at k ^ {1 — 6)w increases 
the amplitude of the tail of the incident wave as we can observe in Fig. HI It contributes so 
relevantly as the peak of the incident wave to the final composition of the transmitted wave 
and creates an ambiguity in the definition of the arrival time. 



IV. ONE DIMENSIONAL POTENTIAL SYMMETRICAL SCATTERING FOR 

E<Vo 

In order to recover the scattered momentum distribution symmetry conditions for ac- 
curately applying the SPM, we assume a symmetrical colliding configuration of two wave 
packets traveling in opposite directions. By considering the same barrier represented in 
(HUj) . we solve the Schroedinger equation for a plane wave component of momentum k for 
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two identical wave packets symmetrically separated from the origin x = 0, at time of collision 
t = —{mL)/{2ko) chosen for mathematical convenience, where we assume they perform a to- 
tally symmetric simultaneous collision with the potential barrier. The wave packet reaching 
the left (right )-side of the barrier is represented by 

/ + 00 
dk g{k - k,)(l)^^''^ {k, x) exp [-i E t] (56) 

where we assume the integral can be naturally extended from the interval [0, w] to the interval 
[— cxD, +oo] as a first approximation. By assuming that 0^(«)(fc, x) are Schroedinger equation 
solutions, at the time t = — (mL)/(2A;o) i. e. when the wave packet peaks simultaneously 
reach the barrier, we can write 



b}'^''\k,x) = eiip[±ikx]+R^''"'\k,L)exp[:^ikx] x < -L/2{x > L/2), 

65^(^)(A;,x) = a^^''\k)exp[^px]+ p^^"\k)exp[±px] - L/2 < x < L/2, 
65'(^)(A;,x) = Tl'^"'\k,L)exp[±ikx] x > L/2 {x < -L/2). 



where the upper(lower) sign is related to the index L{R). By assuming the constraints which 
require the continuity of (p^'" and their derivatives at x = —L/2 and x = L/2, after some 
mathematical manipulations, we can easily obtain 

and 

Th^ik, L) = exp [-^ k L] | exp [p(fe) L] [1 exp [2 z 9(fc, L)]] | 

^ ^ ^ 1 l-exp[2/9(A;)L]exp[ie(/c,L)] J ' ^ ' 

where Q{k,L) is given by the Eq. (l43l) and I^{k,L) and Tg{k,L) as well as R'^{k,L) and 

T^{k, L) are intersecting each other. By analogy with the amplitude addition procedure we 



have adopted in the multiple peak decomposition scattering such a pictorial configura- 
tion obliges us to sum the intersecting amplitude of probabilities before taking their squared 
modulus in order to obtain 

{k,L)+T, (k, L) = exp [-z L] | ^ ^ ^^^^^ ^^^^^^ j = exp L - ^{k, L)]} 

(59) 

with 

r 2kp{k) smh.\p{k)L\ \ 

ifik, L) = - arctan <^ — ^'^ !, , } . (60) 

' ^ + {k^ - p^{k))cos\i[p{k)L\\ ^ ' 

The important information we get from the relation given by (!59|) is that, differently from 
the previous standard tunneling analysis, by adding the intersecting amplitudes at each side 
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of the barrier, we keep the original momentum distribution undistorted since the modulus 
'^(fc, L) +Tg'^(k, L) I is unitary. At this point we recover the most fundamental condition 
for the applicability of the SPM which allows us to accurately find the position of the 
peak of the reconstructed wave packet composed by reflected and transmitted superposing 
components. 

The phase time interpretation can be, in this case, correctly quantified in terms of the 
analysis of the new phase ^{k^ L). By applying the stationary phase condition to the recom- 
posed wave packets, the maximal point of the scattered amplitudes g{k — k(^\i^'^{k,V) + 
Ts'^(A;,L)| are accurately given by k^^^ = k^ so that the traversal/reflection time or, more 



generically, the scattering time, results in 



2m L w'^ sinh (a) -\- ak^ ^^^^ 



/cq a 2kl — + w'^ cosh (a) 

with a previously defined. It can be said metaphorically that the identical particles rep- 
resented by both incident wave packets spend a time of the order of inside the barrier 
before retracing its steps or tunneling. In fact, we cannot differentiate the tunneling from 
the reflecting waves for such a scattering configuration. The point is that we have introduced 
a possibility of improving the efficiency of the SPM in calculating reflecting and tunneling 
phase times by studying a process where the conditions for applying the method are totally 
recovered, i. e. we have demonstrated that the transmitted and reflected interfering ampli- 
tudes results in a unimodular function which just modifies the envelop function g{k — k^) by 
an additional phase. The previously observed incongruities which cause the distortion of the 
momentum distribution g{k — kg) are completely eliminated in this case. At the same time, 
one could argue about the possibility of extending such a result to the standardly estab- 
lished tunneling process. We should assume that in the region inside the potential barrier, 
the reflecting and transmitting amplitudes should be summed before we compute the phase 
changes. Obviously, it would result in the same phase time expression represented by (16TI) . 
In this case, the assumption of there (not) existing interference between the momentum am- 
plitudes of the reflected and transmitted at the discontinuity points x = —L/2 and x = L/2 
is purely arbitrary. Consequently, it is important to reinforce the argument that such a 
possibility of interference leading to different phase time results is strictly related to the 
idea of using (or not) the multiple peak (de) composition in the region where the potential 
barrier is localized. 
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In order to illustrate the difference between the standard tunneling phase time and the 
alternative scattering phase time t!^ we introduce the new parameter n = kl^^^/w'^ and the 

classical traversal time r = (mL) / in order to define the rates 

tr 2 f cosh (a) sinh (a) — Q n (2n — 1) 1 . 
= _ = [4„(i_^)+,i^h^(„)] 1 (62) 

and 

r a [ 2?i — 1 + cosh (q) J 
which are plotted in the Fig.([2]) for some discrete values of n varying from 0.1 to 0.9, from 

which we can obtain the most common limits given by 

lim {RtM] = lim {Rria)} = (64) 

and 

Hm {i?T(a)} = 1 + — , lim \lC(a)] = 1 + - (65) 

Both of them present the same asymptotic behavior which, at first glance, recover the 
theoretical possibility of a swj»er/wmma/ transmission in the sense that, by now, the SPM can 
be correctly applied since the analytical limitations are accurately observed. At this point, it 
is convenient to notice that the superluminal phenomena, observed in the experiments with 
tunneling photons and evanescent electromagnetic waves j^, 0, 37, 38|, has generated a lot 



of discussions on relativistic causality. In fact, superluminal group velocities in connection 
with quantum (and classical) tunnelings were predicted even on the basis of tunneling time 
definitions more general than the simple Wigner's phase-time 3J] (Olkhovsky et ai, for 
instance, discuss a simple way of understanding the problem j3]). In a causal manner, it 
might consist in explaining the superluminal phenomena during tunneling as simply due to 



a reshaping of the pulse, with attenuation, as already attempted (at the classical limit) 39| . 
i. e. the later parts of an incoming pulse are preferentially attenuated, in such a way that 
the outcoming peak appears shifted towards earlier times even if it is nothing but a portion 
of the incident pulse forward tail [5|, |40|. In particular, we do not intend to extend on the 
delicate question whether superluminal group-velocities can sometimes imply superluminal 
signalling, a controversial subject which has been extensively explored in the literature about 



the tunneling effect (j25| and references therein). 

Turning back to the scattering time analysis, we can observe an analogy between our 
results and the results interpreted from the Hartman Effect (HE) analysis jsl . The HE is re- 
lated to the fact that for opaque potential barriers the mean tunneling time does not depend 
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on the barrier width, so that for large barriers the effective tunnehng-velocity can become 
arbitrarily large, where it was found that the tunneling phase-time was independent of the 
barrier width. It seems that the penetration time, needed to cross a portion of a barrier, 
in the case of a very long barrier starts to increase again after the plateau corresponding 
to infinite speed proportionally to the distance ^|. Our phase time dependence on the 
barrier width is similar to that which leads to Hartman interpretation as we can infer from 
Eqs. (16411651) ■ Only when a tends to we have an explicit linear time-dependence on L given 
by 

2 m L / 1 \ , , 

= 1 + - 66 

w \ n ) 

which agree with calculations based on the simple phase-time analysis where It = 
(l + 2^). However, it is important to emphasize that the wave packets for which we 
compute the phase times illustrated in the Fig. ([5]) are not effectively constructed with the 
same momentum distributions. The phase <d{k,L) appears when we treat separately each 
momentum amplitude g{k — kf) |T(fc, L)| and g{k — kf)\R{k^L)\ and the other one (p{k,L) 
appears when we sum the amplitudes on the same side of the barrier. The summed ampli- 
tude corresponds to a recomposed symmetrical momentum distribution for which the SPM 
will give the time-position of the peak of the corresponding wave packet. In this sense, for 
overcoming the abovementioned incongruities, we have reconsidered the use of the multi- 
ple peak decomposition technique previously presented in the study of the above barrier 
diffusion problem jl^ 



V. CONCLUSION 



In this manuscript we have discussed the wave packet multiple peak decomposition in 
order to identify some incongruities due to erroneous application of the method of station- 
ary phase for deriving phase times related to wave packets scattered by a unidimensional 
potential barrier. The main point on which we have elaborated shows us that the results of 
the SPM depend critically upon the manipulation of the amplitude prior to the application 
of the method. The method is inherently ambiguous unless we know, by some other means, 
at least the number of separate peaks involved. By analyzing the above barrier diffusion 
problem, we have demonstrated that the barrier results can be obtained by treating the 
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barrier as a two-step process. This procedure involves multiple reflections at each step and 
predict the existence of multiple (inflnite) outgoing peaks. 

In the same analysis which generically involves analytically-continuous gaussian pulses, 
holomorphic functions which do not have a well-deflned front, we have quantifled some as- 
pects inherent to the tunneling phenomenon. We have observed that the amplitude of the 
transmitted wave is essentially composed by the plane wave components of the front tail of 
the incoming wave packet which reaches the flrst barrier interface before the peak arrival. 
As we have noticed, this filter effect by itself leads to ambiguities in the interpretation of 
the wave packet speed of propagation. To overcome these incompatibilities which leads to 
the misconception of the tunneling phase time, we have considered the multiple peak de- 
composition technique in a suitable way for comprehending the conservation of probabilities 
for a very particular tunneling conflguration. An example for which, we believe, we have 
provided a simple but convincing resolution, since the asymmetry presented in the previous 
case was eliminated, and the phase times could be accurately computed. 

We let for a subsequent analysis the suggestive possibility of investigating the validity of 
our approach when confronting with the intriguing case of multiple opaque barriers 35[, in 
particular, in the case of non-resonant tunneling. Still concerning with the future theoretical 
perspectives, the symmetrical colliding conflguration also offers the possibility of exploring 
some applications involving soliton structures. In summary, the above discussion reinforces 
the assertion that the investigation of wave propagation across a tunnel barrier has always 
been an intriguing subject which is wide open both from a theoretical and an experimental 
point of view. 
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1,0 1,5 2,0 2,5 8,5 9,0 9,5 10,0 

k / w 



FIG. 1: The phase derivative 8'(fc) dependence on k/w. Q'{k) does not present an adequate analytical 
behavior (smoothness) for the applicability of the SPM when k approximates to w since the phase derivative 
oscillates too rapidly (The phase is not stationary) . The method can be accurately applied for larger values 
of k/w when the phase is really stationary. 
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TABLE I: The values of k numerically obtained by increasing the barrier extension 
L.The values are calculated in terms of the wave packet width a for different values 
of the potential barrier high expressed in terms of w a. We have fixed the incoming 
momentum by setting ko a equal to 1 . 



w a 

L/a 




1.5 




2.0 




4.0 




6.0 




8.0 




10 




20 


0.00 


1 


0000 


1 


0000 


1 


0000 


1 


0000 


1 


0000 


1 


0000 


1 


0000 


0.05 


1 


0062 


1 


0188 


1 


1777 


1 


4156 


1 


6238 


1 


7726 


1 


9834 


0.10 


1 


0235 


1 


0648 


1 


3799 


1 


6769 


1 


8547 


1 


9397 


2 


0051 


0.15 


1 


0489 


1 


1223 


1 


5349 


1 


8251 


1 


9505 


1 


9937 


2 


0133 


0.20 


1 


0794 


1 


1825 


1 


6571 


1 


9178 


2 


0000 


2 


0204 


2 


0203 


0.25 


1 


1129 


1 


2420 


1 


7575 


1 


9813 


2 


0317 


2 


0390 


2 


0272 


0.30 


1 


1478 


1 


3001 


1 


8430 


2 


0289 


2 


0562 


2 


0551 


2 


0342 


0.35 


1 


1836 


1 


3565 


1 


9185 


2 


0679 


2 


0779 


2 


0704 


2 


0413 


0.40 


1 


2196 


1 


4116 


1 


9874 


2 


1025 


2 


0986 


2 


0857 


2 


0484 


0.45 


1 


2558 


1 


4657 


2 


0524 


2 


1350 


2 


1191 


2 


1012 


2 


0556 


0.50 


1 


2921 


1 


5194 


2 


1155 


2 


1668 


2 


1399 


2 


1170 


2 


0628 


0.55 


1 


3285 


1 


5729 


2 


1785 


2 


1988 


2 


1611 


2 


1331 


2 


0701 


0.60 


1 


3649 


1 


6266 


2 


2429 


2 


2314 


2 


1828 


2 


1495 


2 


0775 


0.65 


1 


4015 


1 


6809 


2 


3101 


2 


2651 


2 


2051 


2 


1663 


2 


0850 


0.70 


1 


4383 


1 


7360 


2 


3819 


2 


3002 


2 


2281 


2 


1834 


2 


0925 


0.75 


1 


4751 


1 


7920 


2 


4599 


2 


3367 


2 


2518 


2 


2009 


2 


1001 


0.80 






1 


8489 


2 


5466 


2 


3751 


2 


2761 


2 


2188 


2 


1078 


0.85 




* 


1 


9065 


2 


6456 


2 


4154 


2 


3013 


2 


2371 


2 


1155 


0.90 




* 


1 


9646 


2 


7627 


2 


4578 


2 


3272 


2 


2558 


2 


1234 


0.95 




* 




* 


2 


9091 


2 


5028 


2 


3540 


2 


2750 


2 


1313 


1.00 




* 




* 


3 


1137 


2 


5504 


2 


3818 


2 


2947 


2 


1392 



° For the values of L marked with *, we can demonstrate by means of Eqs. (|54II55|) that the 
modulated momentum distribution has already been completely distorted and the maximum 
loses its meaning in the context of applicability of the method of stationary phase. 
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FIG. 2: Erroneous interpretation of tlie seattering of an ineoming wave packet by a unidirnensional potential 
barrier of height Vq and width L. It has the purpose of illustrating the deficiencies inherent to the wrong 
applicability of the SPM. We have plotted the propagating wave packets at the corresponding times tn = 
{ma^)[n {L/a)]/{aqo) (with n = 0, 1, 5 and with the normalization constraint ma^ = 1) by assuming the 
incoming wave packet starts at a; = — {ko L)/{2qo). Prom the false behavior of the density of probabilities it 
becomes obvious that the total probability is not conserved as it was expected. The square of the amplitude 
modulus would supposedly represent a collision of a wave packet of average width a with a potential barrier 
Vo of width L = 5a where, for illustrating reasons, we have adopted ko = \/2w and wa = 10000. Only a 
fixed region in x close to the barrier is shown. 
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FIG. 3: Multiple peak decomposition of a propagating wave packet in the above potential barrier scattering 
problem. We have plotted the first few reflected and transmitted waves for times t„ — {•ma?)[n {L/a)]/ {a Qq) 
(with n = 0, 1, 5 and ma^ = 1) in correspondence with Eq. (|34p . by assuming that the incoming wave 
packet starts at x = — (fco L)/{2qo). The density of probabilities represents the collision of a wave packet of 
average width a with a potential barrier Vq of width L — 5a. Just for illustrating reasons, we have adopted 
ko — (\/T0w)/3 {wa — 10000) and we have printed the wave packet amplification multiplying factor (A) 
(individually adopted for visual convenience for each wave packet) when necessary. 
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FIG. 5: Time rates for tlie standard tunneling and the new scattering process. The rates R{a) and R'^{a) 
can be understood as transmitted times in the unities of the classical propagation time r. Both of them 
present the same asymptotic behavior. 
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